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Abstract 



We investigate holographic superfluids in AdSd+i with = 3, 4 in the non-backreacted approx- 
imation for various masses of the scalar field. In = 3 the phase structure is universal for 
all the masses that we consider: the critical temperature decreases as the superfluid velocity 
increases, and as it is cranked high enough, the order of the phase transition changes from 
second to first. Surprisingly, in = 4 we find that the phase structure is more intricate. For 
sufficiently high mass, there is always a second order phase transition to the normal phase, no 
matter how high the superfluid velocity. For some parameters, as we lower the temperature, 
this transition happens before a flrst order transition to a new superconducting phase. Across 
this flrst order transition, the gap in the transverse conductivity jumps from almost zero to 
about half its maximum value. We also introduce a double scaling limit where we can study 
the phase transitions (semi-)analytically in the large velocity limit. The results corroborate 
and complement our numerical results. In c? = 4, this approach has the virtue of being fully 
analytically tractable. 
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1 Introduction 

High-Tc superconductivitjl^ is a challenging theoretical problem because it involves a strongly 
coupled non-linear phenomenon, while theoretical physics is mostly about solving (weak per- 
turbations of) linear equations. One way we might actually be able to solve high- Tc supercon- 



"'^ Typically high- Tc superconductors are Copper oxides of heavy transition elements, called the Cuprates. 
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ductivity is if we were able to translate the strongly coupled problem into a weakly coupled 
different problem, so that we can get to work on it with the usual hammer and tongs. 

The AdS / CFT correspondence [H [21 E] offers such a translation for many strongly coupled 
systems in terms of a weakly coupled theory that involves gravity in anti-de Sitter space. It 
has recently been found that some of the generic features of superconducting phase transitions 
and condensates can in fact be reproduced using an Einstein-Maxwell theory coupled to a 
charged scalar field in an AdS-black hole background [5|, [S]. This gravitational system has 
instabilities towards the formation of scalar hair, and this has a dual interpretation in terms of 
a superconducting/superfluid phase transition. 

Superconductivity implies that the conductivity is infinite for DC currents, which in turn 
means that the system can sustain a steady DC current in the absence of any applied electric 
field. So any gravitational system that purports to be dual to a superconductor should have 
states which can be identified as the duals of such constant current states. The AdS / CFT map 
claims that the boundary values of bulk-fields are sources for the dual fields on the boundary. 
A current Jx along one of the spatial directions x on the boundary is therefore dual to a gauge 
field Ax turned on in that direction. So the holographic dual of a supercurrent would be an 
AdS black hole with vector hair (which should be thought of as the excited current carrying 
state) on top of the hairy scalar black hole ground state (the superconducting vacuum). 

In [01 [7] such a holographic supercurrent /superfluid solution was constructeclfl in AdS^ when 
the charged scalar field was conformally coupled. The purpose of this paper is to generahze it 
to various scalar masses and also consider AdS^ as well. The mass-squareds of the scalars we 
consider range from positive, down to those that saturate the Breitenlohner-Freedman bound. 

In (i = 3 we find that the system exhibits the same phase structure in all cases: in every case 
that we have investigated, the critical temperature for the superconducting phase transition 
decreases as one increases the superfluid velocity 0, and in fact changes from a second order 
to a first order transition at a "special point" in the phase diagram. This was the behavior 
found in [6l [7] for the conformally coupled case @. We find the same general behavior in all the 
cases where we have done explicit direct numerical simulation. These cases all have m? < 4 
(we set the AdS scale to unity). We also find evidence based on a semi-analytic/semi- numeric 

^The bulk U{\) gauge field is dual to a global [/(I) on the boundary. So these solutions are perhaps more 
accurately described as superfluids, not superconductors. But for many physical questions, the fact that the 
C/(l) is gauged is irrelevant because the effective "dressed" mass of the electrons is large. The flow of current 
in a real superconductor is often nothing but superfluid flow. 

Other related works include the non-abelian case [lU], inclusion of gravity backreaction [IS], fixed super- 
current case [n], rotating cases [M], supercurrent vortices [21] and inhomogeneous case [20] , 

''This is the leading order (constant) piece in the falloff of at the boundary. See Section 2. 

^See [H for an interesting phase diagram in a similar anisotropic system, where a backreacted non-abelian 
holographic superconductor was considered. 
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approach that the same is true for masses as high as = 20 around which the numerics become 
unstable. Our results should be contrasted with the results of fTS] where backreaction effects 
where taken into account in the conformally coupled case in d = 3. It was found that when the 
charge of the scalar is small enough (which is the opposite limit to where one can ignore the 
backreaction effects), the phase transition stays second order even at high superfluid velocity. 
Another complementary study is [TT] where conformally coupled case has been discussed in 
canonical ensemble, i.e. in a fixed current ensemble. 

In d = 4, our results are much more surprising. For low mass the situation is similar to 
d = 3 case. However for sufficiently high mass, we find that there is always a second order phase 
transition to the normal phase, for arbitrarily high superfluid velocity. In some intermediate 
regimes, as we lower the temperature, this transition arises at a higher temperature than a 
subsequent first order transition to a new superconducting phase. We describe this structure 
(which we call the Cave of Winds for reasons to do with the shape of the condensate curve) in 
more detail in the main text. We have checked these claims by direct numerical simulation of 
the condensate (and the free energy) for some values of the masses. We have also developed 
an approach that works in a scaling limit at high velocities, which is fully analytical in d = 4. 
The results using that approach are valid for phase transitions from the normal phase. They 
indicate a second order phase transition at very high velocities (for high enough mass), and are 
consistent with our direct numerical results where they overlap. We can analytically determine 
the critical mass at which this change in phase structure happens, and it turns out to be 
approximately ^ —2.457. In Fig. [H we present schematic diagrams of the various phase 
transition scenarios we have found. 

The next section serves the dual purpose of introducing the problem and describing the 
setup. Section 3, which is the core section of this paper, presents out numerical results as well 
as presents the details of the phase structure. We also devote a part of the section to developing 
and applying the (semi-)analytical approach that can be used to study the phases in a high- 
velocity/low-temperature scaling limit. Section 4 discusses aspects of the conductivity plots. 
An appendix that recapitulates some of the relevant aspects of phase transitions is provided. 



2 Setup 

We will work with the black brane in AdS (dimension arbitrary for now) coupled to a 
Maxwell field and a charged scalar. In general dimensions, the matter piece in the action takes 
the form 



S 



I 




^F^'F,t,-V{\tP\)-W-tqAi^\' 



(2.1) 
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Superfluid velocity 



1st Order 




Temperature 

(a) For d = 3 and d = 4 with lower values of m^. The undot- 
ted curve is the line of second order transition which ends in a 
special/critical point. The dotted line is for first order transition. 



Superfluid velocity 



Normal 




Phase III 



2nd order 



Temperature 

(b) Phase diagram with d = 4 and intermediate values of 
m^. At high enough superfluid velocity, there is a flrst order 
transition between superconducting phase III and super- 
conducting phase I (dotted line). There is always a second 
order transition between normal phase and phase I (undot- 
ted line) . The behaviour shown near T = is tentative. 
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Normal phase 



Phase I 



2nd order 



Temperature 

(c) Phase diagram for = 4 with large values of m^. Only a 
second order transition is possible between superconducting 
and normal phase. 



Figure 1: Phase structure 
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We have set the gauge couphng g = 1, which is merely a choice of unit. But the charge q of 
the scalar is still arbitrary. We would like to take a limit where we can treat the matter part 
of the action merely as a probe, so that the matter piece is negligible compared to the Einstein 
piece in the full action. To do this, we redefine "if = qip and A = qA, so that the matter action 
takes the form 

S = -^J d'^+^x^ (^--^F'^'Fab - Vm) - - z^^p) . (2.2) 

Now it is clear that sending g — )■ oo while holding \E',^ fixed takes us to the probe limit. In 
particular, the matter action dies relative to the Einstein action, so the full equations of motion 
for the metric will be the vacuum equations. We can then solve for ^>,A in the given Einstein 
vacuum solution, because the matter EOMs are unaffected by the overall factor of 1/q^ in the 
action above. By writing the complex scalar as = ip e*^ and doing the gauge transformation 
Ab = Ah — Vf,6', we can get rid of the phase and bring the equations of motion to the form 



V'Fab = 2V>2^, (2.3) 
V^V,^ - 2(VMj7/> - 2^A"V,^ - (A"Aj7/> - ^\/'(^) = 0. (2.4) 

The metric background we will work with is the planar black hole (AdS black brane): 

ds^ = -fde + f~Hr^ + r^dxidxi, f{r) = (^1 - ^) . (2.5) 



We make the choice that the AdS radius is set to unity. The boundary of AdS is at r — )■ oo in 
these coordinates. To turn on a supercurrent, following [HIE] we look for a bulk vector potential 
of the form Aa = {At{r),0, Ax{r),0, ...), where the dots indicate that all of those components 
are zero. The equations of motion in the above background take the form 




We will set V{iIj) = m^il)'^ from now on, and try various values of the mass in what followqj. 
With this choice, the potential term in the EOMs becomes linear. 

For a generic mass m? the falloff of the scalar and the gauge field at the boundary can be 
determined by solving the asymptotic form of the equations above, which are nothing but the 

^In d — Z the choice = — 2 corresponds to the conformally coupled case considered in [51[7]. 
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corresponding equations in pure AdS. One then gets the following falloffs: 

^ = ij^ + i^ + ..., where A± = -(d ± Vd^ + 4:m'^), (2.7) 
r-^- r + 2 

A, = /.--^ + ..., A^ = S^--^^+.... (2.8) 

For the scalar, the usual AdS/CFT correspondence ^ identifies the boundary value of the 
normalizable mode to a vev in the dual field theory and that of the non-normalizable mode to 
a source that couples to the vev (we discuss some subtleties in this statement momentarily). 
Our aim is to set the source term to zero and compute the vev by starting from the horizon and 
integrating out to the boundary. For > —d'^/4+1, only the A+ mode is normalizable, so there 
is no ambiguity: we set = and identify = (O). In the regime —d'^/4 < m? < —d'^/4 + 1 
both modes are normalizable, so we can choose to set either to zero, and treat the other as 
the condensate. In the plots, we report both. For the boundary case = —d'^/A where the 
scalar is at the Breitenlohner-Freedman bound, the two A's coincide and therefore there is a 
logarithmic branch. We set this logarithmic branch to zero because otherwise it triggers an 
instability [9]. 

A+ A_ 

-4 2 2 
3z Z I 

4 2 2 

_ 15 5 3 
4 2 2 

-3 3 1 
"O 4 0~ 

Table 1: The cases we consider in ci = 3 (left) and d = 4 (right). We list the dimensions of both 
the normalizable and non-normalizable modes, but set the non-normalizable mode to zero to 
do the numerics. When both modes are normalizable, we consider both possibilities separately. 

In each dimension, the minimal choices would be to consider (1) the case corresponding to 
the BL bound, (2) a case in the range where both modes are normalizable, and finally (3) a 
case where only one mode is normalizable. The cases we consider in gory detail are collected 
in Table 1. Besides these, we have also looked at some special cases less exhaustively. We also 
introduce a semi- analytic approach in a high velocity scaling limit, and consider a range of 
masses in both d = 3 and d = 4. This has the advantage of being simpler, but at the expense 
of loosing some information regarding the detailed structure. In d = 4 near ^ —2.5 the 
phases have a lot of structure and we have not attempted an exhaustive scan. In both d = 3 



2 

m 


A+ 


A_ 


9 
4 


3 
2 


3 
2 


-2 


2 


1 





3 
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and (i = 4, we consider the conformally coupled scalaiu in detail. In d = 3, this is the case with 
= —2 and was already investigated in [HI E]. 

The problem as we have stated so far has a scaling symmetry that we can use to simplify 
the numerics. We note that we can rescale 

r^r/ro, At At/ro , A^ ^ A^/vq , ^ ip , Xi ^ Xi/ro , (2.9) 

to get rid of the horizon location tq from the equations. In other words, we set the horizon 
location at r = 1. We also rescale the boundary coordinates to keep the boundary metric in 
the usual Minkowski form. 

The gauge field component At has to be zero in order for the norm of A to be finite at the 
horizon [S] . Expanding the fields in a Taylor expansion around the horizon r = 1 and plugging 
into the equations of motion using At{r = 1) = we find the constraints 

^' = ^{Al-m')i;, (2.10) 

K = l^^i^', (2.11) 

at the horizon. 

We are interested in the dimensions of the various quantities in terms of the temperature 
of the black hole. The temperature T goes linearly with rg, so by using ( 12. 9p we find from (12. 7[ 
12.81) the dimensions of the various quantities as follows: 

[^_] = A_, [^+] = A+, [f,] = [S4 = l, [p] = [J,]=d-l. (2.12) 

Since /j, has the dimensions of temperature, T//i is a dimensionless quantity. As we have chosen 
to rescale ro to unity, we will work with l//i instead of temperature. Before we proceed, a word 
about terminology: we call the superfluid velocity and the current. This is reasonable, 
because in conventional condensed matter physics, the gradient of the phase of the condensate 
is what is usually called the superfluid velocity [17]. In turn, the gradient of the phase can 
be traded for the vector potential by a gauge transformation as discussed previously. So it 
is natural to associate the leading piece in the falloff of the vector potential to the superfluid 
velocity. Furthermore, the superfluid velocity can be thought of as the source of the supercurrent 
as dictated by AdS/CFT P]: the Lagrangian for the superfluid [IT] indeed has such a form. 

generic scalar Lagrangian in a curved background is not Weyl invariant. But if the Lagrangian consists 
of the standard kinetic term and a mass-hke term with "mass" of the form £,d+iR where R is the Ricci scalar, 
the action happens to be Weyl invariant if £,d+i = ^^^^jj^- This is called a conformally coupled scalar. When 
the background is AdSd+i, R = — i^i^-tll (we have retained the AdS scale), and the "mass" = — -^^^^r^ is a 
constant and becomes a true mass. When we set L = 1 this gives m? = —2 in d = 3 and to^ = "X in d = 4. 



7 



The simple strategy in what follows is to integrate the EOMs from horizon to boundary. 
These are three second order equations, so we need six numbers to fix them. The horizon value 
of At being zero and the resulting constraints fl2.10p fix three of them. Setting the source term 
in the scalar fixes another. So we are left with a two-parameter family of solutions. This means 
that there are secret relations between (O) = ip+,T = 1/fi, and J^- What we are interested 
in is the plot of (O) in the S'^,T plane and see how robust the phase structure is, as we tune 
the mass of the scalar and the dimension of spacetime. 

3 Phase Structure 

In this section, we present the results of the computations in d = 3, 4 for various values of the 
mass and for various values of the superfluid velocity. As explained in the previous section, we 
will use n to set the scale in the system. This means that we are working in the grand canonical 
ensemble, where the charge density (as measured in fi) can vary, while the chemical potential 
(as measured in fi) is fixed. We can work also in the canonical ensemble, where p is what sets 
the scale of the system in which case we will be in the canonical ensemble. The general nature 
of the plots that we present are qualitatively unaffected by this switch (even though we have 
checked this explicitly only for a handful of cases). Each curve is for fixed S = Sx/p, while 
scanning temperature, T = When presenting the final result, we have found it useful 

to normalize the temperatures in terms of the critical temperature at zero velocity, T^. The 
representative values of S that we have chosen are noted on each plot. It is also important 
to note that the condensate axes are also appropriately scaled (as indicated on the plots) for 
convenience. To be precise, for an operator Ox of dimension A that condenses, we plot {Ox)/T^ 



Following [51 E] we study our system in a grand canonical ensemble, i.e. we keep the 
superfluid velocity fixed when comparing between different solutions. In principle one may 
work in the canonical ensemble and keep the supercurrent fixed [TT] . In the paper by Horowitz 
and Roberts [12], various condensates of holographic superconductors were investigated. We 
have checked that our results precisely agree with their results when the superfluid velocity is 
zero. 

3.1 Instability and the Normal Phase 

The hairless solution to the equations (12. 7|) is given by 



vs. T/T,. 



A: 



■X 




(3.1) 
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This is the normal phase solution that exists for all temperatures. This solution is just the 
non-backreacted version of the RN solution {i.e. the backreaction of the gauge fields on the 
geometry is ignored). Our goal here is to argue that as we increase fi, and thus decrease 
T = the scalar field in the above background develops a tachyonic mode and condenses. 
The result can be thought of as a new bound state in an auxiliary Schrodinger problem as we 
now explain. 

The scalar EOM in (12. 7p may be written as (assuming V{iIj) = m^ip'^) 

.,, f f d—l\,, /A} Al\ , m? , , , 

r + + -fjY - = • (3-2) 

After applying the following change of variables: 

ip dr 1 

r— dy f 

the scalar EoM (13. 2[) takes the form of a Schrodinger equation: 

-^i,-v,siy)i' = o. (3.4) 

Notice that ?/ — )■ oo as r — )■ 1 and ?/ — )■ as r — > oo. The potential Ves, written in terms of r, 
reads 

,f (d-l)id-3) {d-l)f , Al Al m^ 
M0 = -/ ( ^^T^-7rf-—)- ^3.5) 

The EOM of ^l) has been rephrased as a Schrodinger-type potential problem on a semi-infinite 
line {%) : [0, oo]). Depending on the nature of the potential \4fr, there may exist a bound state 
for tl). A bound state such as this signifies an instability, suggesting that there may be new 
phases with non-trivial The resulting new phase (with non-trivial scalar and gauge field) is 
called a "superconducting" phase (and the solution is called a holographic superconductor) as 
it shows infinite DC conductivity [8]. 

In the absence of A^ and At^ the scalar ■?/' may condense if and only if rn? < m\p. This 
is the well known Breitenlohner-Freedman (BP) bound. Generically, by choosing a sufficiently 
negative V^g one may force ip to condense, e.g. if we set 5*^. = 0, = yu(l — i) and increase 
from zero then eventually the solution will become unstable. For even larger values of /i, there 
will be a bound state of ip. The existence of such a bound state may be rigorously shown using 
a trial wave function [16j. The situation is a little more complicated if Sx 7^ 0. It is clear from 
the nature of the potential that a non-zero discourages condensation. Using the method of 
trial wave function it may be possible to address this problem, but we will not attempt that in 
the present paper. 

Conceptually, however, it is fairly easy to understand that a zero mode of tp is only possible 
if and only if S' = ^ < 1. Plugging the solution (13. ip into the eq. (13. 5 p the potential takes the 
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form: 




^) , (3.6) 



Kff(r) = ^^V + V;%, where V = - [{l - - S'^) 



(3.7) 



and Vgjj is the effective potential for the scalar in the absence of gauge fields. Considering the 
fact that V is a simple polynomial in 1/r, it may be shown that for S" > 1, V is non-negative. 
Hence V^g > V^^j and we conclude that the lowest eigenvalue of in the potential V^g is strictly 
greater than the lowest eigenvalue of ip in the potential V^g. In particular, Ves cannot have a 
bound state if V^g does not already have one. Then for 5 > 1 there cannot be an instability 
towards the condensation of ip. In the case < 1, it is easy to see that the lowest value of V 
is S*^ — 1, occurring at the boundary r — )• oo. Moreover, for sufficiently large values of /x, V^s 
becomes negative enough in a region of finite support near the boundary. Hence it is natural 
to expect a condensation of ip. These observations are indeed supported by our numerics. 

Unfortunately, the EOM of is not exactly solvable in general. Numerically one may find 
out the trajectory in the {S, ^) plane where the zero mode of ip occurs. However, even with 
Sx = the EOM of tp is not solvable in general. Although an exact solution near the phase 
transition for = —4 in d = 4 has been presented in [18|, till date there is no known exact 
solution with non-zero Sx- Consequently there is not much analytical understanding of the 
nature of the phase transition at finite Sx- Here we will show that taking a certain double 
scaling limit enables us to exactly solve the EOM for ip. We will thus prove some interesting 
analytic results about the phase transition using this scaling limit. But before doing that we 
present the detailed numerical solution of eq. fl2.7p for various specific cases and discuss the 
corresponding phase structure. 



The phase transitions in the d = 3 case, for all masses that we have investigated, follow the 
pattern that was already observed in the conformally coupled = —2 case [01 E] • At high 
temperature the only possible phase is the normal phase. At small S the superconducting 
phase transition is just like the Sx = case of [8]. The normal phase becomes locally unstable 
at T = Tc (not to be confounded with Tc which all along this paper stands for the critical 
temperature at Sx = 0) and a new superconducting phase (I) is created for T < T^. The 
superconducting phase has less free energy than the normal phase and dominates over it. This 
fact may be analytically argued for Sx = 0- The associated phase transition is of the second 
order. 



3.2 d = 3 
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Figure 2: S = 0.0, 0.32, 0.75 increas- 
ing to the left. 




Figure 3: S = 0.0, 0.36, 0.64 increas- 
ing to the left. 




Figure A: S ^ 0.0, 0.24, 0.40 increas- 
ing to the left. 



j,3/2 




Figure 5:5 = 0.0, 0.30, 0.52 increas- 
ing to the left. 
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d = 3 , = -9/4 , A = 3/2 





d = Z, = -9/4 , A = 3/2 




Figure 6: A typical first order phase transition: S — 0.64 for d — S^m^ — —9/4. On the left we 
plot the condensate versus the temperature and on the right we show the corresponding free 
energy zooming in on the region close to the phase transition. 
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Tc decreases with increasing superfluid velocity and above some critical value of 5* the order 
of the phase transition changes from second to first. This is reflected in the free energy plot 
(see Figure [6])§| as the transition from a smooth curve to a swallowtail cusp§|. As we lower 
the temperature, before Tc is reached 0, two superconducting phases (apart from the normal 
phase) become possible at a temperature Tjq. Unlike the previous case both these phases have 
a non-zero below T/y. The phase with larger value of ip (phase I) has less free energy than the 
phase with lower value of ip (phase II). The free energy of both superconducting phases (I and 
II) decreases as T is lowered and at some temperature, say T = Ti, the free energy of phase I 
becomes lower than the one of the normal phase. Consequently phase I becomes the dominant 
phase in the system and continues to be so for T < Ti. The associated phase transition is of 
first order. At a smaller temperature T = T^. (note that < Ti) the phase II merges with the 
normal phase. 

In a later section, we show by a semi-analytic scaling limit approach that the high-velocity 
first order transition is present for masses as high as = 20. We have done detailed direct 
numerical simulation for a handful of masses (m^ < 4), some of which we present here in the 
condensate curves and the free energy plots of figures [21 [31 IH E] and [HI The phase structure is 
quite robust. 



3.3 d = 4 

The d = 4 case is more interesting. Here, we find that as the mass of the scalar increases, 
qualitatively new features appear. In particular, for high enough mass, the phase transition 
remains second order for arbitrarily high values of , as is evident for example from our = 
plots in figures [S] and [71 For small enough mass the phase structure is similar to that of the 
d = 3 case. In the intermediate mass scale the phase diagram shows interesting new features. 

"Cave of Winds" : A numerical scan of all values of the mass is quite challenging, so we 
have settled for plotting the graphs for a few representative values. A particularly interesting 
representative case is = —7/4 (see figure [Hj). Here, for low values of 5* the phase transition 
is always second order. For high enough superfluid velocity, we find that as the temperature is 
lowered, there is at first a second order phase transition from the normal to a superconducting 
phase and then a first order phase transition between two superconducting phases (see figures 
[TOl and [HD. We call this the Cave of Winds phase structure because of the shape of the 
condensate plot. The basic structure here as we lower the temperature is that initially two new 

^In this plot and the ones of figures [Tj and [TT] we have subtracted the free energy of the normal phase. 

discussion of the free energy is presented in the section [5^ 
^"In the cases with many "critical" temperatures, we use Tc to denote the phase where the normal phase 
becomes locally unstable. 
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phases (phase I and phase II) appear, both of which are of higher free energy than the normal 
phase. As we further lower the temperature, a new superconducting phase (phase III) with 
lower free energy branches out from the normal phase. At this stage, phase III is the dominant 
phase. The transition from normal to phase III is second order. At a still lower temperature, 
there is a first order transition and phase I becomes the dominant phase. It is remarkable that 
the overhang of the condensate curve stretches farther than the location where the transition 
happens: this is is the reason why the first order transition is to another superconducting phase. 
A phase transition that was further along the overhang would result in a first order transition 
directly to the normal phase, and would be in contradiction with the analytical scaling limit 
results that we report in a later section. Happily, this is not the case. A further comment that 
is worth making regarding the relation between the analytic scaling limit and direct numerical 
simulation is as follows. The scaling limit works in a high velocity, low temperature limit and 
it sees the phase structure close to the normal phase. So it is not surprising that it misses 
the first order transition between the two superconducting phases that we see via the direct 
numerical approach. It sees only the second order transition. Another comment is that (as we 
observed above) the first order phase transition happens long before the tip of the overhang: 
we have noticed that this is a fairly generic phenomenon that happens for other cases as well, 
and not just for the Cave of Winds case. In the language of Ginzburg-Landau effective actions, 
a cave of wind like situation could appear if some higher order term (higher than quartic) in 
the effective action changes sign: we discuss this in the appendix. 

For the even higher mass case of = 0, the phase transition is always second order for 
any superfluid velocity. In a later section we investigate various masses in c? = 4 using a high 
velocity scaling limit, analytically. The results suggest that the critical mass at which the phase 
transition structure changes is ^ —2.457. Below that the phase transition at high velocities 
is first order. The analytical approach suggests that above that mass the higher temperature 
phase transition is second order at high velocities, but it is silent about what happens at lower 
temperatures. Indeed, for the values of mass that we have done detailed numerical simulation, 
we find that at high velocities there is always a second order transition. While the analytical 
scaling limit is rather simple and can be applied to many cases at once, it suffers from the 
drawback that it cannot be applied to understand phase transitions away from the scaling 
limit. This means in particular that it only sees phase transitions in which at least one of 
the phases is the normal phase. Therefore it has nothing to say about (say) the first order 
transition that we see in the = —7/4 case. Of course, in all the cases where we expect 
overlap between the two approaches, we find matches. There seems to be a lot of structure in 
the intermediate range of masses close to ^ —2.457 and we hope to come back to a detailed 
numerical investigation of this range in future work. 

Finally, in figures [121 1131 HH and [15] we show the results of our numerical analysis for the 
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To 



Figure 7: A typical second order phase transition: S — 0.97 for d = A^m^ — 0. On the left 
panel we plot the condensate versus the temperature while on the right panel we present the 
corresponding free energy, zooming in on the region close to the phase transition. 




Figure 8: S ^ 0.0, 0.61, 0.8 increasing Figure 9: S ^ 0.0, 0.8, 0.9 increasing 

to the left. to the left. 

cases = —3, —4 and —15/4 in d = 4. As mentioned above, for these low values of the mass 
we find the same phase structure as in the d — 3 case. 

3.4 Free Energy 

To understand the nature of the phase transitions, we now describe the computation of the free 
energy. We will plot the free energy of the normal phase and the superfiuid phase for various 
values of the superfiuid velocity to corroborate our claim that the phase transition changes 
from second order to first order as we increase S^- 

The free energy of a solution is defined essentially as the on-shell action after the appro- 
priate addition of boundary (counter) terms: Q = —TSqs- Substituting our ansatz in terms of 
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d = 4 , = -7/4 , A = 7/2 




0.130 0.135 0.140 0.145 0.150 0.155 0.160 0.165 T 



Figure 10: Zoom on the S = 0.9 case of figured The vertical (red) hne indicates the first order 
phase transition between two superconducting states characteristic of the Cave of Winds phase 
structure. 
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Figure 11: Free energy (and its zoom) for 5* = 0.9, c? = 4, w? = —7/4. 
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Figure 12: S = 0.0, 0.36, 0.74 increas- 
ing to the left. 



Figure 13: S = 0.0, 0.33, 0.55 increas- 
ing to the left. 
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Figure 14: S = 0.0, 0.32, 0.68 increas- Figure 15: S = 0.0, 0.35, 0.56 increas- 

ing to the left. ing to the left. 

At{r), Ax{r),ilj{r) and the background (12.51) in the action of the Einstein-Maxwell-scalar theory 
given by eq. (12. 2p we get 

^bare = - [ d'^'x T^-' f-4^ + ^-,^1 + - tA'^' + + , (3.8) 



J V 2 2r2 ^ / r2 ^ 

and after using the EoM's (12.71) we arrive at 

(3.9) 

Notice that we have both boundary and bulk terms contributing to the on-shell action and the 
bulk contributions are always finite in the cases we consider here. From the falloffs (12. 7[ 12.81) . 
we find that the only divergences possible are from the scalar terms in the boundary action. 
The scalar terms are schematically of the form 

where r is the radius at which we cut-off the geometry!"]. The first piece is divergent and is 
killed off by the counter-term [TTl [T3]. The final piece is zero when we take the cut-off to 
infinity. The finite piece is the middle piece and its coefficient gets corrected due to the finite 
part of the counter-term. Even the remaining finite part will not contribute to the free energy 
in our case because one of the two among ip^ is always zero for us. The specific counter-term 
is chosen on the basis of the boundary conditions and ensemble that we are working with, see 
[7] and section 3.1.1 of [H] for a discussion. In our case here, we are working with superfluid 
velocity fixed as opposed to the current fixed. 



^^This expression is written for the case when the A's are not equal. But the resuh for the final expression of 
the finite part holds even in that case. 
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The upshot of the discussion above is that the scalar pieces do not contribute in our com- 
putation of the free energy. Then, the final expression of the free energy that we will use in our 
numerical analysis is 

d r 



47r 

1 («) A («) 



d-2 



(Af Af-'^ - 4-^-2)) - dr( r^-'ij'Al - 



f 



,(3.11) 



where {Ai , Ai ) correspond to the falloffs of the dimensionless rescaled gauge fields 



A, 



id-2) 



+ ... 



A - 



(d-i) 



rd~2 



+ ... 



(3.12) 



and are related to the dimensionful falloffs of eqs. fl2.7[l2.8p through the rescaling (see eq. (12. 9p ^ 



/i = ro Af^ , p = ry Ar~'' , = vo , J. = r^' A^'^ . (3.13) 



d-l A{d~2) 



The T"^ and the overall factor {An/dY in fl33T|l arise precisely because we are working with 
dimensionless variables to perform the numerics. We have also divided by the volume of the 
boundary. As we have mentioned in footnote 8 in the free energy graphics of figures El [7] and [TT] 
we plot AQ: the free energy of the phase with condensate minus the free energy of the phase 
without condensate, i.e. the normal phase corresponding to the analytic solution of eq. (13.11) . 



3.5 Nature of the phase transition 

Let us try to have some analytic understanding of the nature of phase transition. From eq. 
(13.81) by integration by parts and using the scalar EoM (13. 2p we get 



n 

J'd 



+ ; 


[ d'^+'x ( 


r=oo J 





. f A 12 



d'^+^x r^-^ + ) , (3.14) 



2 2r2' 



where in the last line we have used the fact that the boundary value of is kept to zero in our 
case. Writing the action in this simple form gives us information about the relative free energy 
of the different phases. 

The free energy for the normal phase is given by '^'y™"' = —{d — 2)^/i^/2, as one can easily 
check by plugging the hairless solution (13. ip into eq. (I3.14p . As discussed before, as we increase 
the value of /i, condenses and there is a new superconducting phase. The outstanding 
questions are: (1) Which phase dominates? (2) What is the nature of the associated phase 
transition? For Sx = ^ this question may be answered without getting into too much detail. 
With a chemical potential fixed to /i, the solution of At in the new phase can be written as 

A =/^(l-^) + 5A, (3.15) 
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where 6 At — )■ at r — )■ 1 and r — )• oo. Then, from eq. f l3.14p we get, 



W 2 

,2m 0^2 



+ {d-2) j drdr{5At) - j drr" 

oo f 

+ {d- 2)5 At ^ - h 



12 

t) 



^^V-2Y , ,,,, - fdrM^ 



'2 



2 

(3.16) 



Hence 

= - f^now = -{T'^V) j dr r"-^^!^ < . (3.17) 

Therefore if a phase with non-trivial ip exists it will always have a lower free energy than the 
normal phase. 

A little extension of the above argument determines the nature of the phase transition 
locally. As we increase /z for fixed a zero mode of ■?/' = ■i/'o forms for /i = /ic. Following the 
philosophy of [l8l|19], we may consider turning on a small amount of zero mode and at the next 
order look at the backreaction of this mode on At and A^. This enables us to calculate the free 
energy of the new phase with non-trivial At and A^. Let us start by choosing ip = eip^ + o(e^), 
where e is a small constant. This perturbation takes the system away from the normal phase 
to the superconducting phase. We can plug in this value of ip and solve for 5 At and 5Ax with 
appropriate boundary conditions. It is clear from the equations that the variation of At and 
Ax will be proportional to e^, i.e. 

At = Af^ + ^5At + o{e'), 4°^=^(l-^. 
Ax = Sx + eHAx + o(e^) 

where 



(3.18) 



dri 



SA. = S.j j^-^j dr2^\'-\ (3.20) 

We fix the constants of integrations by the boundary conditions at the horizon, i.e. 6 At 
vanishes and SA^ is regular. It should be noted that boundary values of At and also changes 
by o(e^). With respect to the normal phase, the relative action of the present configuration is 
given by 



W 2 



j drr'-'(^{d6Atf-^{d6Ax)')+{d-2)6fi'] + o{e') 



Se^ + o(e^) . (3.21) 
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The sign of S determines tlie pliase transition. If 5 < 0, tlie associated transition is locally 
second order. If = 0, this is the only possibility as we saw before. On the contrary, the 
situation becomes more interesting if 5 > 0. This happens for a typical first order transition. 
Although our small fluctuation analysis does not address the question of global phase structure, 
a positive S strongly suggests (from a simple Ginzburg -Landau picture) the existence of a first 
order transition (see appendix). Notice that we have some more information we can use: for 
fixed Sx, a very large value of fi implies that the effect of would be neghgible and we expect 
a similar phase diagram to the A^ = case. Therefore the dominant phase in the large /i limit 
is the superconducting phase with non-trivial ip. But since, as we have just explained, the 
superconducting phase never becomes locally dominant over the normal phase, the only way in 
which it can become dominant is through a first order transition. 



3.6 Double scaling limit and exact solutions 

In this section we will show that in a suitably defined double scaling limit where — ?■ 1 and 
/i — oo, the EOM of ip simplifies drastically. In this limit the potential Ves defined in eq. (13. 5p 
is negative only in a very small neighbourhood near the boundary y = (or r = oo), so all the 
interesting condensation mechanisms occur there. Therefore we now expand eq. (13.41) around 
y = and keep only the leading terms: 

^ ^ ^ + /i^ ((1 - 2y'-') -S')^ + {m' + ^^^)y'-'i, = , (3.22) 

where we have taken into account that at leading order in /(y) ~ p-, = /i(l — y'^"^), and 
= 'S'a;. Let us now define the following double scaling limit: 

S — > 1 , /i — > oo , 

il = /i2-2°(l - S^) kept fixed and a = (3.23) 

a 

In this limit equation fl3.22p reduces to 

— ^ + (/i - 2z'~y = 0, (3.24) 

where we have introduced z = yfi"". This equation happens to be exactly solvable in c? = 4 for 
general mass, while the conformal case (m^ = — -^^-j^) is also solvable in ci = 3. We shall first 
focus on the d = 4 case. 



^^Note that the first correction to / is of o{y'^) and subleading with respect to the y dependent part of At 
so it can be ignored consistently. On the other hand, the y dependent part of At is necessary to generate a 
non-trivial potential for -ip and survives our double scaling limit. 
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3.6.1 d = 4: 



In d = 4, we can solve eq. fl3.24p exactly. Generically this equation may be solved in terms of a 
product of Laguerre polynomials and Gaussians. However we will not present the full solution 
here. Our interest is confined to finding a node-less, normalizable solution. Such a solution is 
given by 

4,{z) = z~^exp{-^), (3.25) 

with ^ 

A(A-l) = m2 + ^-^. (3.26) 

In addition, we must require the following relation between A and fl: 

fl^ = V2{2\ + 1). (3.27) 

It might seem surprising that we are demanding that there be a constraint on the parameters 
of the differential equation that we started with. The appropriate comparison here is with 
the energy E in the standard Schrodinger equation of quantum mechanics, which can get 
quantized due to boundary conditions. It should be noted that eq. f l3.24p does not always have 
normalizable solutions. These only occur for the specific values of fi{= fie) given by the above 
relation. These are then the values of fl for which we may find zero modes of ip. 

As follows from eq. f l3.26p A has two solutions for each value of m^: 



A = A± = -(1 ± V4m2 + rf2). (3.28) 



It should be kept in mind that as follows from eq. (13. 3 P ^jj ~ y^i^ near the boundary (y ~ 0). 
Consequently, using eq. (I3.25P t/^ ~ z'^ ~ z^~^^ (at z ~ 0). Notice that A = A + (ci — l)/2 is 
the dimension of the operator dual to and has the two possible values written in eq. (12. 7p . 
The unitarity bound gives the constraint A > — | |22]. Thus A+ is always a valid choice for A 
and for — ^ < m? < — ^ + 1, A_ is also a valid choice. In the tabled we have summarized the 
properties of the system for a few specific values of A. 

It is straightforward to address the eigenvalue problem associated with eq. (I3.24p . For 
jl > fic, the zero mode discussed here becomes a negative energy bound state and ip tends to 
condense. No such instability exists for /t < /ic. 

Besides the ground state there exist excited solutions with nodes in the z direction. These 
solutions are given in terms of Laguerre polynomials. We will not explore them in full detail 
and only write down the solution for = — 

1 ^2 

^(z) = H„(24z)exp( -), f={2n + l)V2, n e Z* . (3.29) 

Furthermore, A = 1 implies an odd n and A = implies an even n. 
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3.6.2 d = 3 



Here we will look at the conformal case i.e., m? = —2. Normalizability at 2; = oo fixes 

^ = AiryAi(^^^), (3.30) 

up to a multiplicative constant. There may be two possible boundary conditions near the 
boundary z = 0. We may choose either ip{0) = or ip'i^) = 0. They correspond to the 
condensation of the two possible dual operators. Choosing the appropriate boundary condition 
fixes the value of ^ for which a zero mode of ip exists. 

The boundary condition ip{0) = implies 

fic = 2ho, (3.31) 

where zq is some zero of the AiryAi(z). It should be noted that all zeros of AiryAi(2;) lie in the 
2; < region. Choosing the first zero (nearest to 2; = 0) gives us a solution which does not have 
a node. Choosing other zeros give rise to multi-nodal solutions [23]. It is believed that these 
solutions always have higher free energy than the node-less soliton and do not play a significant 
role in phase transition. Choosing the first zero of AiryAi(2) we get 

^(0) = ^ /ie ^ 3.71151 (3.32) 

On the other hand, the boundary condition ijj'{0) = implies 

/ic = 2i2o, (3.33) 

where Zq is some zero of the AiryAiPrime(z) = ^AiryAi(2;). Similarly to the previous case, 
all zeros of AiryAiPrime(2;) lie in the 2; < region. The first zero (the one nearest to 2; = 0) 
gives us a solution which does not have a node. Choosing other zeros gives rise to multinodal 
solutions as before. For a node-less solution 

^'(0) =0^/1^^ 1.61723 (3.34) 



3.6.3 Free energy in the double scaling limit 

Here we will concentrate mainly on the d = 4 case. We will follow the general philosophy of 
section 13.41 to calculate the free energy in a small fluctuation analysis. Our calculation would 
be in the double scaling limit (13.231) where S ^ 1. In the scaling limit we may simply take 
= 1 (we are in the region near the boundary) . Changing to the radial variable y and 
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keeping the leading order terms in e we get from eq. fl3.2ip 
5n = V^eVi / dyy'^-^ 



dx ip'^ 



id -2) 



dyy' 



d-3 



dx - x''-^) 



(3.35) 



which in terms of the rescaled variable z = yfi^ becomes 



2V e n d < He dz z 



dx ip 



2 dz z' 



d-3 



dx x'^ "^ip"^ / dyip 



(C/-2) 



dz z 



d-3 



dx ip' 



(3.36) 



Here we have kept only the leading terms in fi. This expression is exact in the double scaling 
hmit 0323]). 

In c? = 4 the above integral may be evaluated analytically using the exact solution fl3.25p . 
We get 

2 



sn 



-3A- 



4M4P-l)r( A + l ) -8Ar(2A + l) 



(3.37) 



The Gamma function can be expressed in terms of integers, radicals and vr's when A is an 
integer or a half integer. In table [2] we summarize a few values of A and the corresponding 
values of 5Vt. We also plot the dependence of 5VL on A numerically (figure [T6|l . For small values 
of A, 5VL is positive and it becomes negative for large values of A. It may also be argued using 
Sterling's approximation that 5VL is a negative quantity in the large A limit. 5VL changes sign 
between A = | and A = 2, to be precise at A ~ 1.742 (corresponding to rn? ^ —2.457). 

As we discussed before, a positive value of 6Q indicates that locally the superconducting 
phase has more free energy than the normal phase and there is no second order transition from 
the normal to the superconducting phase. However this opens up the possibility of a global 
phase transition between the normal and superconducting phases. Considering that for fl > fic 
the normal phase is locally unstable, there possibly exists a global transition before ft reaches 
flc- On the other hand, if SQ is negative then locally the superconducting phase has lower 
free energy than the normal phase and there will be a second order phase transition from the 
normal to the superconducting phase. 

The situation is more complicated in the d = 3 case. Only a piece of the integral (I3.36P 
can be computed exactly. In particular, 6fi may be calculated exactly, but the other second 
integrals can only be reduced to first integrals. Unfortunately an analytic evaluation of the 
resulting first integrals seems unlikely and the final result has to be calculated numerically. 
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sn 




Figure 16: Dependence oi 6Q on X. A negative value of 6Q implies a second order transition 
between the superconducting and the normal phase. As we discussed before a positive SQ gives 
a strong hint for a first order transition in the system. 
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Table 2: In this table we display the values of the free energy for several condensates in c? = 4 
in the double scaling limit. In the first two columns we present the dimension of the operator 
and the mass of its dual bulk field, while in the third column we show the corresponding value 
of the free energy computed using eq. fl3.37p . The free energy changes sign between A = 3 and 
A = |. 
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4 Transverse AC Conductivity 



In this section we discuss the AC conductivity results for two characteristic cases: one where 
the phase transition changes to first order at high velocity and another where it remains second 
order all the way. We will only consider the transverse conductivity here: the longitudinal 
conductivity couples to the other modes and things are more complicated. The idea is to solve 
the transverse gauge field perturbation 6 Ay = e^^'^^Ay{r) numerically in the background with 
the condensate. The equation for it takes the form 

~ 2 d— 3 

dr{r^-^fdrAy) + 2r''-''i,^)Ay = 0. (4.1) 



We have defined 



where the last step uses the relationship between the Hawking temperature of the black brane 
and its Schwarzschild radius. According to linear response theory, the conductivity is defined 
by a Kubo formula in terms of retarded current-current correlators, which via the AdS/CFT 
dictionary can be related to the boundary falloffs of Ay 

•4. = + ^ + - (4.3) 
The final result for the conductivity [12] in c? = 3 is 



A 



(1) 
y 



uA^' 



and in d = 4 



2A(2) 



The above result (and the retarded correlator) is obtained via the Schwinger-Keldysh pre- 
scription by solving the equation with infalling boundary conditions at the horizon. We do this 
numerically and plot the results for two representative cases in figures [TTJ and [181 The first 
corresponds to the case d = 4, m? = (A = 4), which is a case that exhibits second order phase 
transitions at all velocities. We plot the real part of the conductivity versus the frequency 
for different values of the superfiuid velocity at a fixed temperature T = 0.22Tc. At large 
frequencies the behavior is linear while for small enough frequencies Re(cr) vanishes within our 
numerical precision, this is the conductivity gap. As expected, when the value of the superfiuid 
velocity is close to its maximum (where the phase transition happens) the gap vanishes. This 
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Re(cr 




Figure 17: Plot of Re((7) for d = 
4, = at temperature T = 

0.22Tc. The velocities are S = 
0.000278, 0.427, 0.708, 0.840, 0.850, in- 
creasing (on the lower part of the fig- 
ure) to the left. The highest veloc- 
ity corresponds to the point where the 
2nd order phase transition to the nor- 
mal phase occurs. 




Figure 18: Plot of Re(cr) for d = 
4, = —15/4 at temperature T = 
0.16Tc. The velocities are S = 
0.000320, 0.505, 0.759, 0.905, increas- 
ing (on the lower part of the figure) to 
the left. At the highest velocity the 
system undergoes a 1st order phase 
transition to the normal phase. 



is consistent with the second order nature of the phase transition, and should be contrasted to 
the second case (figure [TS]) where we plot the conductivity for d = A^im? = —15/4 (A = 5/2 
and T fixed to 0.16Tc). Here, as we increase the velocity the gap decreases, until we reach the 
maximum velocity where the gap still has a finite value, consistent with the first order nature 
of the phase transition at high velocities. Roughly, a first order phase transition means that 
the condensate has a finite binding energy while second order transition corresponds to the 
case where the condensate can have arbitrarily low binding energy: so the nature of the phase 
transitions we observe is consistent with the results for the frequency gap. 

The study of the AC conductivity will lead us to an interesting observation about the Cave 
of Winds phase structure we discussed in section 13.31 Let us indeed consider the Cave of Winds 
condensate plots and take two points on either side of the first order phase transition. Note 
that both these are superconducting states: the one on the upper branch corresponds to what 
we called phase I before and the one on the lower branch corresponds to phase III. The results 
for the conductivity are presented in Fig. ( 1T9|) where we clearly see that during the first order 
phase transition the conductivity gap falls from half of its maximum value (corresponding to 
T — > 0) to almost zero. It would be very interesting to see if there are any real material that 
exhibits such a smoking gun for a "Cave of Winds" structure: as one lowers the temperature 
in the superconducting phase, the conductivity gap would jump from almost zero to half the 
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2 4 6 S ^ 

Figure 19: Plot of Re((j) for d = 4,m^ = —7/4 at fixed velocity S = 0.9. The different lines 
correspond to T equal to (from left to right): 0.147Tc (phase III), 0.147Tc (phase I) and 0.102Tc 
(phase I). 



value at zero temperature. 



5 Summary and Comments 

Our original motivation for the work was that the construction of [6l[7] was undertaken in the 
probe approximation where the charge of the scalar q is large and the backreation on the black 
hole geometry was ignored. This is a somewhat more drastic assumption in the supercurrent 
case than in the case of the superconducting ground state: the perturbation breaks some of the 
symmetry of the background, and to take account of a fully backreacted solution we will need 
to consider a non-static (but stationary) metric. So it is interesting to see whether the results 
of O E] where robust. 

During the course of our work, [15] appeared, which solved the backreaction problem for 
the conformally coupled case in = 3, and it was found that for small enough values of the 
charge the phase transition stays second order for all velocities (at least to within numerical 
control). In the non-backreacted situation that we are considering, we find another kind of 
deviation from the results of [6l [7] : when we go to = 4, we find that the phase structure can 
take quite an interesting structure for large enough mass. We saw the emergence of a sequence 
of phase transitions (of different order) with temperature and also the possibility of a new 
superconducting phase. Another interesting direction would be to look at zero temperature 
limit [23]. 

It would be interesting to investigate the phase diagram in more detail in the region of 
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masses near = —2.457 in = 4, where there seems to be a lot of structure. It might be 
possible that, beyond what we have already seen, extra features show up in these cases. A plot 
of our phase diagram where the mass is an extra axis would be interesting because it is not 
clear how the Cave of Winds phase diagram transitions into the first order structure that we 
see at high velocities at lower masses. This would require a rather exhaustive numerical scan 
of the phases. 

Another obvious line of development is to consider the backreacted version of some of these 
cases where one can ask low-temperature questions with more control. This would basically 
involve the construction of various kinds of fully backreacted hairy black holes in AdS^. In 
flat 3+1 dimensional spacetime, hairy black holes are ruled out by the famous results of Israel, 
Carter and Robinson, as well as general arguments by Hawking. In 4+1 dimensions on the 
other hand, by now it is well known that there is continuous hair In AdS as well, it seems 
likely from our results that the situation in d = 4 could be richer than what is possible in d = 3. 
Yet another natural question is that of the embedding in string theory. Finally, it would also 
be very interesting to see if there are any real systems in condensed matter which show the 
kind of Cave of Winds structure that we saw. In particular, the jump in the gap will be an 
interesting feature of such materials. The fact that we saw this in d = 4 (and not in d = 3) 
suggests that these will probably have to be 3+1 dimensional systems. We hope to come back 
to some of these problems in the future. 
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6 Appendix 

A. Effective action and phase transition 

Here we review some elementary aspects of phase transitions. Let us consider a simple generic 
effective action 

S = ^<P' + l<P' + l<P'^ OO, (A.l) 
with an order parameter 0. We choose c to be positive to prevent a runaway situation. 
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Let us start with a positive and gradually lower it. For a positive m^, the trivial solution 
ip = would be locally stable. As changes sign the trivial solution becomes locally unstable. 
The mere fact that the trivial solution becomes locally unstable does not automatically imply 
a second order transition. The order of the transition depends on the sign of 6. If 6 > the 
phase transition is locally second order. However if 6 < 0, then a first order transition happens 
before changes sign. This feature does not strictly depend on the particular form of our 
effective action and only relies on the fact that the effective action has no runaway directions 
and is not unbounded from below. 

At the phase transition point there is a zero mode of ip. In our simple case the zero mode 
is 00 = 1- Near the phase transition point we may turn on a little amount of the zero mode 
around the trivial solution to perturbatively construct a new solution. Let us assume that 
= e + o(e^). Then from the EoM of we find 



Hence the on-shell action of the new phase is related to the sign of b and consequently determines 
the nature of phase transition. If 6 < (or 6Q > 0), the perturbative phase never becomes 
a local minimum. Hence when the trivial phase becomes locally unstable, the minimum of 
free energy lies in some other minimum well separated from the trivial phase. (In the case of 
a runaway situation such a minimum exists for — oo.) This explains why the associated 
transition would be more abrupt than a second order transition and would generically be first 
order. 

Also to be noted is the parameter range at which the new perturbative phase exists. For 
6 > the new phase exists for < 0, while for 6 < it arises for > 0. 

To rephrase the above argument in the language of our problem, we look at the off-shell^ 
effective action of the lowest eigenvalue mode of ip, schematically written as. 



It may be argued that for small /i, is a positive number. At /i = /ic, changes sign. 
Generically near = fic, = i^{fic — /^)- If 5 < a new dominant phase is created for 
fi > ^c- However if 5 > there is no new phase for fi > fic- 

^^Note that it is the zero mode of the auxiUary Schrodinger problem that corresponds to the on-shell case. 
Here we are considering a non-zero eigenvalue. 





Hence the on-shell action for the new phase is. 




(A.3) 




(A.4) 
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The Cave of Winds Structure : 



A Cave-of- Winds-like structure is possible if we include the effect of higher order terms in our 
effective action. Let us consider an effective action of type 

S-'^^' + l^' + l^' + i^"^ ^>0 (A.5) 

Now if c changes sign and becomes sufficiently negative, while b remains positive, then there 
will be a Cave of Wind like structure. Locally near cj) — the phase transition will always be 
second order. However because c is negative, there may exist a globally dominant phase which 
becomes dominant through a first order transition. Whether this phase transition happens 
before or after the second order transition {i.e. when crosses zero) is a question of detail. 
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